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O ' ABSTRACT 

CO ! We use the joint measurement of geometry and growth from anisotropic galaxy clustering in 

IT) . the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) CMASS sample 

' reported bv lReid et alj to constrain dark energy (DE) properties and possible deviations from 

the General Relativity (GR). Assuming GR and taking a prior on the linear matter power spec- 
{NJ ■ trum at high redshift from the cosmic microwave background (CMB), anisotropic clustering 

of the CMASS DR9 galaxies alone constrains Q m = 0.308 ± 0.022 and 100Q k = 5.9 ± 4.8 for 
w = -1, or w = -0.91 +0.12 for Q. k = 0. When combined with the full CMB likelihood, the 
addition of the anisotropic clustering measurements to the spherically-averaged BAO location 
\ increases the constraining power on dark energy by a factor of 4 in a fiat CDM cosmology 

^ • with constant dark energy equation of state w (giving w = -0.87 ± 0.05). This impressive 

gain depends on our measurement of both the growth of structure and Alcock-Paczynski ef- 
fect, and is not realised when marginalising over the amplitude of redshift space distortions. 
Combining with both the CMB and Supernovae Type la (SNela), we find Q m = 0.28 1 + 0.014 
and 1000Q k = -9.2 + 5.0 for w = -1, or w Q = -1.13 + 0.12 and w a = 0.65 + 0.36 assuming 
Qt = 0. Finally, when a ACDM background expansion is assumed, the combination of our 
estimate of the growth rate with previous growth measurements provides tight constraints on 
the parameters describing possible deviations from GR giving y = 0.64 + 0.05. For one pa- 
rameter extensions of the flat ACDM model, we find a ~ 2<r preference either for w > - 1 or 
slower growth than in GR. However, the data is fully consistent with the concordance model, 
and the evidence for these additional parameters is weaker than 2cr. 

Key words: gravitation - cosmological parameters — dark energy — dark m atter — distance 
scale — large-scale structure of Universe 
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1 INTRODUCTION 

All currently available cosmological observations, including the 
latest datasets of CM B temperature and polarisation anisotropies 
dKomatsu et alj 1201 lh . Supernovae Type la (SNela) magnitudes 
l lSuzuki et al.l I2012T) and the distance ladder mapped by Baryon 
Acoustic Oscillation (BAO) p eak signature in the clustering of 
galaxies (Ander son et al.ll2012l) are consistent with a simple cos- 
mological model in which general relativity (GR) describes grav- 
itational interactions on all scales and times, about 70 per cent of 
the Universe's current energy density is in form of a Cosmological 
Constant as originally described by Einstein, and most of the re- 
maining 30 per cent is in fo rm of nonrelativistic "dark matter" (For 
a detailed review see e.g., IPeebles & Ratrall2003l ; Weinberg et al. 
l2012h . 

Ongoing and future observations have been designed to test 
if the cosmological constant needs to be replaced by a dynamical 
Dark Energy (DE), and if so, to measure the properties of this DE 
fluid. We should also be able to tell if GR describes the properties 
of gravity on cosmological scales or if it must be re placed by a yet 
unknown modified theo r y of gravity (MG) (see e.g., Albrech t et al.l 
120091 : IZhao et alj[20Qg| : Isamushia et alj|201ll ; IWang et alj|2010h . 
Observational effects of dynamic DE and MG are partially degen- 
erate an d careful data analysis should take into account both possi- 
bilities dlshak, Upadhve & Spergelll2006l : IShapiro et al.ll2010l) . 

The clustering of galaxies provides a very powerful and ro- 
bust test of both the nature of DE and MG. The shape of the mea- 
1 1 —1 

sured correlation f unction (Reid et al. 2010) or the power-spectrum 
(Montesano . Sanchez & Phlepsll201 1 J), analogously to the shape of 
the measured CMB power- spectrum dLarson et al.l 1201 lb . can be 
used to constrain basic cosmological parameters. Features within 
the clustering signal, particularly the BAO, allow the clustering 
to be used as a standard ruler. Additionally, although the statis- 
tical properties of galaxy clustering are expected to be isotropic, 
the measured clustering can be highly anisotropic, depending on 
how redshifts are translated to distances. The two main sou rces of 
this a pparent anisotropy are redshift-sp ace distortions (RSD; Kaisei 
1 1987b and the Alcock-Paczynski (AP; lAlcock & Paczvnskilll979l) 
effect. 

RSD arise because peculiar velocities contribute to observed 
galaxies redshifts, but can not be corrected for when computing 
line-of-sight separations. On quasi-linear scales, the average pair- 
wise galaxy peculiar velocity is negative, meaning that galaxies are 
on average falling towards the mass over-densities traced by neigh- 
bouring galaxies. These coherent motions appear as a "squashing" 
of the correlation function along the line-of-sight. The amplitude 
of the observed anisotropy can be used to infer the strength of the 
gravitational interaction at differ ent scales a nd redshifts. (For a de- 
tailed review of RSD see Hamilton! fl993) . RSD allow measure- 
ments of the amplitude of fluctuations in the velocity field, which 
in linear theory give a dependence on 



f(z)cr s (z) : 



(1) 



din a 

where cr 8 (z) is the overall normalisation of the matter density fluc- 
tuations. 

The AP effect results from the fact that, to convert observed 
angular positions and redshifts of galaxies into physical positions, 
we must use a cosmological model on the observed lightcone. If 
the wrong model is used when computing the correlation function, 
the initial isotropy of the clustering signal will be distorted. The 
measured anisotropy of clustering can be used to infer the proper 
geometry and hence the true values of cosmological parameters. If 



we have a prior on the shape of the correlation function, the dila- 
tion of scales between the spherically averaged observed and model 
clustering statistics allows a measurement of 



D v (z) = 



1/3 



(2) 



where D A (z) is the physical angular diameter dist ance and H(z) is 
the Hubble expansion rate (Eisenstein et al. 2005). Applying the 
AP test to the measured direction dependent clustering additionally 
allows the measurement of 

1+Z 



m = 



-D A (z)H(z). 



(3) 



This allows the degeneracy between radial and angular distances 
in BAO measurements to be broken and an accurate estimate of 
Hubble expansion rate at different redshifts to be derived. 

Many RSD measurements have been made from a vari- 
ety of galaxy surv eys, includin g most recently the 2dFGRS 
Percival et al1l2004h. t he VVDS fcuzzo et alj|2007h, the 2SLAQ 



da Angela etalj|2008l). the SDSS-II JCabre & Gaztanag j 120091: 
Song et al J201 lbHSamushia et alj2012j), the WiggleZ dBlake et al 
2011al) .andthe6dFGRS ( Beutler et alj|2012l) surveys. These mea- 



surements have in turn been used to set constraints on the cosmo- 
logical growth rate. 

Thus far geometric constraints from galaxy clustering have 
focused predominantly on spherically averaged power spectra or 
correlation functions. From such measurements, the BAO feature 
allows few percent-level distance measurements (for t he latest con- 
l201(t iBlake et alj|201 lbl; iBeutler et al.l 
l2012r . lAndersonetal.1120121) . By corT 



straints, see 



Percival et al 



l201ll ; Ipadmanab han et al 



trast, the AP test has received less attention, since better signal-to- 
noise data is needed to disentangle RSD and AP effects, and more 
careful modelling of anisotropic correlation function is required. 
The AP effect has been recently used to jointly m easure D A , H and 
/erg in three redshift bins from Wiggl eZ survey ( Blake et al .120121) 
and using the SDSS-II LRG sample Ichuang & Wang|l2012l). ~ 

The Sloan Digital Sky Survey (SPSS; lYork et al.ll2000h has 
mapped over o ne third of the s ky using the dedicated 2.5-m 
Sloan te lescope jGunn et alj|2006l) . A drift-scanning mosaic CCD 
camera l|Gunn et ak 19981) i maged sky in five photometric band- 
passes dFukugita et alJI 19961) to a limiting magnitude of r - 22.5. 
The ongoing Baryon Oscillation Spectroscopic Survey (BOSS; 
ISchlegel, White & Eisensteinll2009l) . which is a part of SDSS-III 
jEisenstein et al J|20 1 lh . is measuring spectra of 1.5 million galax- 
ies and 150,000 quasars selected from the multicolor SDSS imag- 
ing. The CMASS sample in the BOSS Data Release 9 (DR9) con- 
tains a catalog of 264 283 highly biased (b ~ 2) galaxies sampling 
an effective volume of about 0.77 h~ 3 Gpc 3 with mean redshift of 
z ~ 0.57, allowing for best-to-date st atistical uncertaint y in the 
measurement of galaxy clustering (see Iwhite et al.lfeoi ll for de- 
tailed description of early data). 

This work is part of series of papers providing a comprehen- 
sive description of the galaxy clustering in the CMASS DR9 sam- 
ple. iNuza et alj J2012h compared the clustering of CMASS DR9 
galaxies to state of the art dark matter simulations a nd showed 
that th ey are compatible with ACDM model and GR. iRoss et al] 
|2012f) provided an up-to-date description of the CMASS DR9 data, 
a study of possible observational s ystematic effects and the meth- 
ods to remove known systematics. iManera et al. I J2012T) presented 
600 mock catalogs that match the observed volume of CMASS 
DR9 sample and are essential for determining the covariance ma- 
trix ass ociated with the mea sured correlation functions and power 
spectra. [Anderson et al.ld2012l) used these data to measure the BAO 
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peak position to 1.7 per cent precision. ISanchez et all d2012l) de- 
rived cosmological constraints using the full sha pe of the measured 
spherically-averaged correlation function, while iReid et all I2OI2I) 
studied the anisotropic clustering of CMASS DR9 galaxies using 
the measured monopole and quadrupole moments of the correlation 
function (henceforth we shall simply refer to these as the monopole 
and quadrupole). 

By adopting a sophisticated mod el for galaxy clustering in the 
quasi-linear regime. IReid et al. U20ll made accurate RSD and AP 
measurements from the direction-dependent clustering of CMASS 
DR9 galaxies to simultaneously measure growth and geometry at 
a redshift of z = 0.57. We now extend this work to investigate the 
cosmological implications of these measurements. We show that 
information provided by the RSD-derived growth rate significantly 
enhances constraints on basic cosmological parameters compared 
to the case where only geometric measurements are used. We com- 
bine measurements of the growth rate, angular diameter distance 
and expansion rate with previous measurements to constrain prop- 
erties of DE and gravity. In combination with CMB, Ho and SNela 
data we are able to estimate values of basic cosmological parame- 
ters to very high precision and tightly constrain possible deviations 
from ACDM and GR (henceforth we will refer to the model in 
which the background expansion follows ACDM and the gravity is 
described by GR as ACDMGR). 

The paper is organised as follows: in Sec. |2]we describe the 
CMASS DR9 AP and RSD measurements, in Sec.[3]we briefly re- 
view previous similar measurements, in Sec. |4]we describe differ- 
ent ways of looking for deviations from ACDMGR, and in Sec. [5] 
we present constraints on deviations form ACDMGR. We conclude 
and discuss our results in Sec. [6] 



2 CMASS DR9 MEASUREMENTS 



Reid ct al. (2012) used the measured anisotropic correlation 
function of galaxies in the CMASS DR9 sample to jointly esti- 
mate D v (z = 0.57), F(z = 0.57) and fcr s (z = 0.57). To obtain 
these estimates they fitted the monopole and quadrupole in the 
range of scales 25Mpc hT 1 < s < 160Mpc hT 1 . The theoretical 
model corr elation functio n was computed from the "streaming 
model" of IReid & White] ( 1201 ll) . Linear theory predictions for 
the anisotropic correlation function cease to be accurate on 
relatively large scales (for the importance of nonlinear corrections 
see and different ways of modeling RSD i n quasi-nonlinear 
regim e see, e.g., Ide la Torre & Guzzol 1201 2t iMatsubaral 1201 it 



Taruva, Nishimichi & Saito 



20 id: 



Taruva, Saito & Nishimichi 



201ltlOkamura. Taruva & M atsubara 201 ll) . IReid & White! ( l20lTh 



demonstrated that the "streaming model" used as the theoretical 
basis of our CMASS DR9 analysis reproduces the monopole and 
quadrupole measured in Nbody simulations with better than 2 per 
cent precision down to the scales of 25Mpc fT 1 ; bett er than the 
level d emanded by DR9 data. Appendices B3 and B4 of lReid et alj 
d2012h further quantify the (small) theoretical errors, including 
uncertainties in the "Fingers-of-God", using realistic mock galaxy 
catalogs. The shape of the linear real-space correlation function 
was marginalised over using a WMAP7 prior on parameters Si m h 2 , 
Q. b h 2 and n s . This triplet of parameters fixes relative clustering at 
different physical separations; the geometrical measurements of 
D\j and F result from matching the measured shape of the galaxy 
correlation function to the shape fixed by CMB data while most 
of information about fa% comes from the r elative amplitudes 
of the measured monopole and quadrupole. IReid et all d2012h 



considered a number of models with different assumptions about 
how the shape of the correlation function relates to the background 
geometry and growth. Most relevant to this paper is the model in 
which the shape of the correlation function is fixed by a WMAP7 
prior but the expansion of the Universe at lower redshifts and the 
growth of st ructure since rec ombination are allowed to be arbitrary 
(Model 4 in lReid et al]|2012t) . Joint estimates of fcr t , D A and H 
recovered by applying this model are independent of assumptions 
about either the late-time behaviour of DE or the structure growth 
and can be used to constrain models with arbitrary DE Equation 
of State (EoS) and growth history as long as physics at recombi- 
nation is unaltered an d linear growth remains scale-independent. 
The IReid et al.l ( 120121) mea surements have already been used to 
constrain Gal ileon gravity ( App leby & Linden k o 1 2 ) and cosmic 
growth index dHudson & Turnbulll2012tlRapetti et al.ll2012h . 

The most general way of testing DE and MG models against 
anisotropic clustering measurements made with the CMASS sam- 
ple would be to fit predictions of the models to the measurements of 
monopole and quadrupole directly. We have checked that in all rel- 
evant cases simply assuming a 3-dimensional multi-variable Gaus- 
sian likelihood around the /<x g , £> v and F measurements of Model 
4 retains most of the constraining power of the data set and re- 
sults in very similar constraints compared with exploring the full 
likelihood. In fac t, while the likelihood of the measurements of 
IReid etal.l ( f2012l) is close to a multivariate Gaussian, it is slightly 
skewed towards higher values of F and fcrg. We have checked that 
the effect of the skewness on the derived likelihood of cosmologi- 
cal parameters is small and only affects regions of low likelihood 
(Appendix [A). For the rest of the paper we will ignore this small 
effect. 

Although th e measurements of geometry and growth rate in 
lReidetal.l ( f2012l) were derived using a WMAP "shape" prior, the 
resulting estimates are not correlated with the WMAP parameters 
Q. m h 2 , Sikh 2 and n 8 . This allows us to combine these measurements 
with the WMAP likelihood without double-counting the CMB in- 
formation. We have checked that combining Model 4 measure- 
ments with WMAP data gives very similar likelihood surfaces as 
those calculated by jointly fitting the CMASS DR9 monopole and 
quadrupole and the CMB data. This approximation allows us to 
greatly speed up the calculation of likelihoods for a variety of cos- 
mological models. 



3 PREVIOUS MEASUREMENTS OF GROWTH AND 
GEOMETRY FROM ANISOTROPIC CLUSTERING 

The growth rate of structure has now been measured from the RSD 
signal in the galaxy clustering pattern, at a range of redshifts from 
z ~ 0.1 up to z ~ 0.7. Fig.Q]shows the most recent CMASS DR9 
estimates of /<x 8 alongside previous measurements from other sur- 
veys that are used later on in our analysis. Measurements, l<x er- 
rorbars and references to original publications are presented in Ta- 
ble[i] 

The bands on Fig.QJshow the oretical pre dictions of GR, PGP 
dDvali, Gabadadze & Porattil200d) and f(R) dBuchdahllfl97cl) the- 
ories of gravity when a WMAP7 prior is applied to cosmological 
parameters describing expansion of the Universe. For the DGP pre- 
dictions we ignore the scale dependence of growth r ate on large 
scales while for the f(R) we use the model studied in Star obinskvl 
d2007l) ; the general f(R) model has a clear GR limit and by tuning 
model parameters its predictions can be made arbitrarily close to 
the GR predictions. 
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0.7 



0.6 



quasistatic f(R) (vertical) 
GR (checkered) 
flat DGP (horizontal) 




Figure 1. The data points show the CMASS DR9 measurement of /erg 
(circle) along with similar, low redshift, measurements (squares) and lcr 
eiTorbars as presented in Table [T] The three stripes show theoretical pre- 
dictions for different gravity models allowing for uncertainty in the back- 

[ round cosmol ogical parameters, constrained using only the WMAP 7 data 
Komatsuetall201ll) . 



fos 



survey 



reference 



0.17 

0.22 

0.25 

0.37 

0.41 

0.57 

0.6 

0.78 



0.51 ±0.06 
0.42 ± 0.07 
0.35 ± 0.06 
0.46 ± 0.04 
0.45 ± 0.04 
0.43 ± 0.03 
0.43 ± 0.04 
0.38 ± 0.04 



2dFGRS 
WiggleZ 
SDSS LRG 
SDSS LRG 
WiggleZ 
BOSS CMASS 
WiggleZ 
WiggleZ 



Percival et al. (2004) 
Blake et al. (2011a) 
Samushia et al. (2012) 
Samushia et al. (2012) 
Blake et al. (2011a) 
Reidetal. (2012) 
Blake et al. (2011a) 



Blake et al. (2011a) 



0.067 
0.77 



0.42 ± 0.06 6dFGRS 
0.49 ±0.18 VVDS 



Beutler et al. (2012) 
Guzzo et al. (2007) 



Table 1. Compilati on of RSD based fcr & measurements. Note that the 6dF- 
GRS measurement Beut ler et al was too recent to be included in our 

analysis, while the measurement of iGuzzo et aTl {2007) was excluded due 
to more recent, stronger, measurements at similar redshift s. 



Geometrical measurements are also available from the local 
Universe up to redshift of z ~ 0.8. Fig. [2] shows the most recent 
CMASS DR9 estimates of D v and F along with similar results 
from other surveys. The numerical values of these estimates, lcr 
errorbars and references to original publications are presented in 
Tables[2]and[3] 

The bands on Fig. [2] show theoretical predictions of the 
spatially-flat ACDM model when a WMAP 7 prior is applied to the 
relevant cosmological parameters, and the Einstein-DeSitter (EdS) 
model (fl m = 1) with Hq = 73.8 ± 2.4km/Mpc/.s (as measured by 
iRiess et aljfoTTh . 



z 


D v Mpc 3 


survey 


reference 


0.106 


456 ± 27 


6dFGRS 


Beutler et al. (2011) 


U.J J 


1 ion -4- 6S 

1 JoU It DO 


jLJjj 1 . I\ V J 




0.57 


2071 ±44 


BOSS CMASS 


Reid et al. (2012) 


0.60 


2234 ± 115 


WiggleZ 


Blake et al. (2011b) 


0.35 
0.57 


1356 ±25 
2094 ± 34 


SDSS LRG 
BOSS CMASS 


Padmanabhan et al. (2012) 
Anderson et al. (2012) 



Table 2. Constraints on Dy from recent surveys shown in F ig. [2] We 
do not include the B AO measurements of Padmanab han et all j2012l) or 
lAnderson et al.1 J2012I) in our analysi s, as they are correlated with the mea- 
surements of Reid et al. (2010, 2012). 



survey 



reference 



0.22 
0.41 
0.57 
0.60 
0.78 



0.28 ± 0.04 
0.44 ± 0.07 
0.67 ± 0.026 
0.68 ± 0.06 
0.97 ±0.12 



WiggleZ 
WiggleZ 
BOSS CMASS 
WiggleZ 
WiggleZ 



Blake et al. (2011c) 
Blake et al. (2011c) 
Reid et al. (2012) 
Blake et al. (2011c) 
Blake et al. (2011c) 



Table 3. AP constraints on F from recent surveys shown in Fig. [2] 



4 MODELLING DEVIATIONS FROM ACDMGR 

The standard cosmological model comprises three main assump- 
tions: first, the Universe on large scales is homogeneous and 
isotropic, secondly, DE is Cosmological Constant and thirdly, grav- 
ity is described by GR on all scales. In this model the background 
geometry of the Universe can be fully described by three numbers 
that can be chosen to be the current relative energy densities of 
nonrelativistic matter Q. m and cosmological constant C1 A and the 
current value of Hubble expansion rate H . The angular diameter 
distance and Hubble expansion rate at any redshift are given by 



D A (Z) : 



dz' 



flb(i + zrUa 

H(z) = H E(z), 



(4) 
(5) 



where 



the quantity fi k = 1 
spatial curvature and 



x if n k = 

— == sinh I Vfik*) if fiw > ,, x 

V"k V ; (O) 

1 sin ( yJ-£hx\ if Q k < 

— O m - fl A is the relative energy density of 



E(z) = V^d + zP + fikd+^ + fiA 



(7) 



When the three basic cosmological parameters are fixed the 
growth of matter overdensities in GR to linear order in overdensi- 
ties follows the growth equation 



d 2 G I dlnH 
+ 2 + 



din a 2 



dG 



din a I din a 



-n m (a)G, 



(8) 



where a = 1/(1 + z) is the scale factor and the growth factor 
G{a) = 6(a)/8(a ini ) shows by how much the overdensities have 
grown compared to some arbitrary initial time a ini . In a spatially- 
flat ACDMGR model the growth rate / = d In G/d In a can be we ll 
approximated by a fitting formula /(a) = fl m (a) ' 55 jPeebles! 1980). 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



z z 

Figure 2. The data points show recent measurements of Dy (left panel, squares) and F (right panel, squares) described in Table 1 along with latest mea- 
surements from CMASS DR9 data (circle) with lcr errorbars. Red stripes show theoretical prediction of spatially-flat ACDM within the uncertainty in basic 
cosmological parameters as m easured by W MAP 7 data; Blue stripes show theoretical predictions of Einstein-deSitter model within the uncertainty in Hq as 
measured by HST key project iRiessetal.l201ll) . 



4.1 Parametrizing deviations from GR 

A large number of viable alternatives to GR have been suggested 
within the scientific community, each with its own theoretical jus- 
tifications and adv antages (for a recent review of MG see, e.g., 
IClifton et alj |2012|) . The deviations from GR are usually con- 
strained in terms of phenomenological parameterisations, with two 
main approaches, either parametrizing the observables (such as / 
and G) or the perturbation equations (Eq.[8). 

The most widely used parameteris ation for the growth rate f 
allows for deviations from GR through dWang & Steinh ardt 1998) 

/ = n m (fl) y , (9) 

where 7 is a free parameter that mildly depends on background 
geometry. In GR this dependence can be approximated well by 

7 = 0.55 + 0.05[1 + w(z = 1)] if w > -1, (10) 
7 = 0.55 + 0.02[1 + w(z = l)]ifw<-l. (11) 

dLindedl2005T) . More accurate expressions for 7 as a function of 
cosmological parameters and redshift that take into account higher 
order terms in 1 - £l m and curvature an d are correct at sub- 
percent level have be en proposed (see e.g., Ilinder & Cahnll2007l ; 
llshak & Do ssett 2009). In GR we expect to find a value of 7 ^ 0.55. 
Higher values of 7 would be indicative of a slower growth of struc- 
ture than in GR and vice versa. 

An alternative way of parameterising deviations from GR is 
to phenomenologically modify the linear equations for dark matter 
perturbations instead of parametrizing the solutions as in Eq. l[9j. To 
describe the most general deviations from GR two scalar functions 
are necessary to modify the temporal and spatial parts of perturbed 
Einstein's equations and many different choi ces for these func- 
tions have been proposed in literature (see e.g., iBertschin get feoOot 



Caldwell, Coorav & Melchiorri 2007; lAmendola, Kunz & S apone 
l2008l ; lAmin. Wagoner & Blandfordll2008l) . 

We consider the parameterisation in which the Eq. ((8) is mod- 
ified to 

d 2 G I dlnH\ dG 3 

77"^ + 2 + -77— IT- = ^ m («)GM«), (12) 
diner y a ma J a ma 2 

where fi(a) parametrizes the time-dependent deviation from GR 
and could also be a function of scale. Eq. d 12b is valid when the 
dark matter is cold and DE does not cluster on the scales of inter- 
est. The function fi(a) is expected to be equal to unity in GR at all 
times and, inspired by the DGP model, the time dependence of fi(a) 
is sometimes parametrized by 

yt/(a) = 1 +fi s a s . (13) 

Large values of fi s correspond to a stronger tim e dependence for 
the effective Gra vitational Constant (see e.g.. ISong et al.l [20 1 1 al ; 
Izhao et aT]|201ll) . Parameter s describes the time-evolution of the 
modification. The s = case corresponds to a simple rescaling of 
the Newton Constant with fi = (G c ff - Ggr)/Gqr while for s > 0, 
the deviation from GR is time dependent and becomes more pro- 
nounced at late times. Larger values of s correspond to increasingly 
recent deviations from GR; therefore fi s is less well constrained for 
larger values of s. 

Besides modifications to Eq. l[8}, MG theories can also in- 
duce a gravitational slip r\ = O/ 1 ?, which is a ratio of longitudi- 
nal to Newtonian scalar potentials Q in the absence o f anisotropic 
sourc es, r\ is expected to be equal to unity (see, e.g., iDaniel et al.l 
2009). The growth rate data alone are not sensitive to gravitational 

1 In Newtonian gauge <I> = 5g|j/2 and *P = iSgJ/2, where 6g^ is a perturba- 
tion to the metric tensor. 
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where 

W(z) = (1 + z) 3 <'+»o+^ exp (j^fj • (16) 

The change to the DE model modifies both the background 
expansion and the growth of perturbations with respect to the 
ACDMGR model. This allows us to use both geometric and growth 
rate measuremen ts to constraint th e model parameters. To illustrate 
the sensitivity of lReid et alj d2012l) measurements of geometry and 
growth to DE model parameters we plot them on fig. |4]along with 
theoretical predictions for different values of w . To derive theoreti- 
cal predictions we set Q.^, and w a to zero and fixed the distance to the 
CMB and amplitude of fluctuations at the last scattering surface. In 
reality each line on fig.|4]corresponds to a band due to uncertainties 
in remaining cosmological parameters such as Q. m h 2 . When these 
uncertainties are properly taken into account the absolute error in 
wo resulting from Dy(z = 0.57)/r s is of order of 0.24 instead of 
0.1 as a naive inspection of fig. [4] may suggest. Fig. [4] also shows 
that changing the value of wq has an opposite effect on /cr 8 and 
F. Since the actual measurements are strongly positively correlated 
their combination results in a much stronger constraint on wq than 
what their individual marginalised errors imply. 



Figure 3. The data points show some of the measurements presented in 
Table[T]with lcr errorbars; the circle denotes the most recent CMASS DR9 
measurement. The thick solid line shows the GR prediction, the thin solid 
lines show predictions for y = 0.45, 0.50, 0.60 and 0.65 (top to bottom), 
and the thin dashed lines show predictions for /jj = 0.2, 0.1, -0.1 and -0.2 
(top to bottom). n m (z = 0) and erg at recombination are kept fixed to their 
WMAP7 best-fit values. 

slip and would have to be combined with an observable such as 
weak lensing observations, that depend on *F, to yield constraints on 
both jUj and t) fu nctions simultaneously (see e.g., ISong et alj20lT3 : 
IZhao et alj2o"Tlh . 

If the amplitude of density perturbations is known accu- 
rately at some earlier time (for example from CMB temperature 
anisotropy data at recombination), and the background expansion 
history is relatively well known, the growth rate measurements at 
low redshifts can be used to strongly constrain deviations from GR. 
Even small modifications to the strength of gravitational interac- 
tion can build up and result in significantly different predictions for 
f(z)crg(z) at late times; the constraints are stronger for models in 
which modifications set in earlier. 

Fig. [3] shows /(z)c 8 (z) predictions for different MG models. 
To help to visualise the constraining power of current RSD mea- 
surements some of the measurements presented in Table [T]are also 
plotted. 

4.2 Parametrizing deviations from A 

Deviations from a Cosmological Constant are most often 
parametrized by w and w/ a , where the time dependence of the EoS 
of DE is approximated as 

Pde z 

=w + w a - (14) 

Pde 1 + z 

JChevallier & Polarskil l200lL Linden 1 120031) where w = -1 and 
w a = give the ACDM limit. In this parameterisation, Eq. ((7} is 
modified to 

E(z) = + z) 3 + "kd + z) 2 + a A W(z), (15) 



4.3 Constraining the DE potential 

A more general way of describing deviations from A is to con- 
strain an effective Lagrangian of the field responsible for the 
accelerated expansion (see e. g., iJimenez. Talavera & Verdell2012t 
iBloomfield & Flanaganll2012T) . If we assume that DE can be effec- 
tively described by a scalar field that is minimally coupled to grav- 
itation, the Cosmological Constant will correspond to a scalar field 
with flat potential V((f>) = const. Observational data can be used to 
constrain higher o rder terms in a Taylo r expansion around the flat 
potential (see e.g. Jjimenez et al.ll2012h . 

Instead of con sidering a general potential V(<f>) as in 
Jimen ez et al, I j2012l) we look at specific scalar field dark energy 
models. We study two representative theories: a slowly-r olling 
scalar field that has a tracker solution (Pee bles & Ratralll988l) and 
a pseu do Nambu-Goldstone boson (PNGB) in an oscillating po- 
tential l IFrieman et al.1 1995). Both of these models can result in 
a late-time accelerated expansion of the Univers e and have been 
shown to be compatib le with current data (see e.g.. lSamushial20o"9l 
iDutta & Sorb o 2008)0 

If the scalar field is minimally coupled to gravity its evolu- 
tion is governed by a Klein-Gordon type equation in the expanding 
Universe, 

. dV(d>) 

cf> + 3H(z)<P - — p = 0, (17) 
dtp 

where V((f>) is the self-interacting potential and the time-dependent 
energy density of the scalar field, 

Pt = \fi + V{<f), (18) 

contributes to the Friedmann equations. 

For the slowly-rolling tracker field, following iRatra & Peebles! 

(1988), we will consider an inverse power-law self-interacting po- 
tential of the form V(0) oc and will assume that it does not clus- 
ter on the scales of relevance. The quantity a is a positive constant 
that parametrizes the steepness of scalar field potential. In the limit 

For a recent review of scalar field DE models see e.e., lLinderl (2008). 
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Figure 4. The data points show the estimates of f(z)o~&(z), Dy/r s and F from iReid et al with lcr errorbars. Solid lines correspond to theoretical 

predictions for different values of ivn when w ;i and £\ are set to zero while the distance to the last scattering surface and erg at recombination are kept fixed to 
their wMAP7 best-fit values. 



of g — > 0, the potential becomes flat and the model approaches a 
cosmological constant; higher values of a correspond to stronger 
time-evolution. The a parameter has been constrained before using 
compilation of d ata sets and fou nd to be less than 0.7 at the lcr 
confidence level ( Samushia 2009J). 

For the PNGB, following iFrieman et al.l dl995l) . we will con- 
sider an oscillating potential of the form V = M 4 [1 + cos(0//)], 
where M and / are free parameters related to the energy scales 
of explicit and spontaneous symmetry breaking respectively. This 
model possesses technical naturalness and can naturally incorpo- 
rate both the small value of cosmological constant and the small 
mass of the scalar field today. Theoretical p redictions of this model 
have been compared to data (see e.g.. IWaga & Friemanl feOOCt 
iNg & Wiltshire1l200ll : iKawasaki. Moroi & Takahashill200ll) ; to be 
consistent with current observations the PNGB field has to be sit- 
ting on the top of oscillatory potential or rolling down very slowly 
and the mass of the scalar field must be larger than ~ M p i/3 
jDutta & Sorboll2008l) . 

Even though the dynamics of scalar field models can often 
be very well described by a pair of wo, w a parameters, fitting data 
directly to physical models has several advantages. The Wo-w a pa- 
rameterisation implicitly imposes a theoretical prior that all values 
of those parameters are equally likely. By directly examining spe- 
cific DE models we can include the information that some values 
of wo and w a are easier to procure than the others. This approach 
also sometimes reduces the number of independent parameters, im- 
proving the constraining power of the data. 



The acceleration of the kinematic expansion of the Universe 
can be described by the "deceleration parameter" q(z), which is 
defined as 

act da(z) 1 + z 



q(z) : 



(19) 



t dz d(z) 

Negative values of 151 indicate an accelerated expansion and vice 
versa. 

The Om parameter, defined as 

(l+z) 3 -l (l+z) 3 -l 

is expected to be constant (equal to the present day 
matter density fl m ) in spatiall y-flat ACDM Universe 
(Sahni, Shafieloo & Starobinskv 2008D and can also be used 
to test acceleration even if the value of £2 m is not known. In a 
universe dominated by DE, O m is expected to be less than the 
critical value irrespective of DE model. 

Finally, the fact that DE becomes relevant very suddenly at 
low redshifts can be tested by looking at the ratio of energy densi- 
ties of DE to nonrelativistic matter 

Pde(z) = H(zflHl-a m {l+z? 

p m {z) n m (i+z) 3 

In Sec. |5.4| we will use the geometric part of our measurements 
in combination with previous similar measurements tabulated in 
Tables [2] and[3]to reconstruct these model-independent quantities. 



4.4 Model-independent properties of DE 



5 MEASURING DEVIATIONS FROM A AND GR 



The ACDM paradigm leads to a number of interesting conclusions 
about the late time Universe: 

• The expansion of the Universe must be accelerating. 

• The Universe at present must be dominated by DE; the contri- 
bution of nonrelativistic matter being subdominant. 

• The Cosmological Constant must be a relatively recent phe- 
nomenon, quickly declining in importance at higher redshifts. 

Whether or not Universe exhibits these properties can be 
checked with observational data without referring to a specific DE 
model. 



We use the corr elated estimates of /tr 8 , D\/r s and F from 
IReid et al.l d2012h to constrain parameters describing deviations 
from AC DMGR. We us e the Monte Carlo Markov Chain (MCMC; 
see e.g., iTannej 1 19961) method to explore likelihood surfaces. 
We do this by making use of publicly available COSMOMC code 
dLewis & Bridiell2002l) that uses CA MB to compute power-spectra 
for th e CMB and matter fluctuations jLewis. Challinor & Lasenbvl 
2000). The COSMO MC software was modifie d to support a time- 
dependent DE EoS (Fang, Hu&Lewisll2008l) . the SNLS dataset 
dConlevetal]|2qilh. the measurements of growth and geometry 
from lReidetal.U2012h and general measurements of growth rate. 
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Figure 5. Confidence levels (lcr and 2tr) for joint fits to parameters y and 
n m . Solid lines show results for CMASS DR9 estimate of growth rate 
only, while dashed lines show results with previous estimates of growth 
rate added. In both cases the growth rate measurements are combined with 
CMB and SNela data. The dotted line shows the expected value in GR. 



5.1 Constraining deviations from GR 

For modified gravity we measure parameters y, or yu s and s, as de- 
scribed in Section l4~n We assume that the background evolution of 
the Universe follows the predictions of spatially-flat ACDM with 
basic cosmological parameters additionally constrai ned by SNela 
dConlevetall201lh and CMB dKomatsu et aDl201lh data. We also 
combine the measurements of iReid et al.l d2012h with previous 
measurements of the growth factor. These measurements are some- 
what heterogeneous. Most were made by as s uming a fi xed fiducial 
cosmo logical model, while Sa mushia et alj J2012h and IReid et al.l 
d2012h marginalize over the shape. Models with varying levels 
of complexity have been used to address the problem of nonlin- 
ear contamination from structure growth and real to redshift space 
mapping. In this section we will ignore these differences and will 
treat the reported measurements and errorbars as a Bayesian likeli- 
hood for fcrg at the redshift of interest. 

Fig.[5]shows two-dimensional constraints on parameters y and 
Q m . Using only the CMASS DR9 measurement of growth rate we 
obtain y = 0.75 ± 0.09 when Q. m and H are marginalised over; 
when previous measurements of the growth rate from Table [T]are 
added we find y = 0.64 ± 0.05, an 8 percent measurement. The 
addition of previous measurements of the growth rate brings the 
best fit value closer to the GR expecte d value of y = 0.55. I n both 
cases GR is about 2<r from the best fit. [Hudson &Turnbulll ( f2012l) 
were able to get a slightly stronger measurement of y by com- 
bining RSD data wi th a low redshift peculi a r velocity data from 
iTurnbull et all d2012h and lDavis et all 3201 lh . iRapetti et alj d2012h 
find a constraint on y at a similar level of precision, although their 
best-fit value is closer to the GR prediction. 

A similar two-dimensional likelihood contour plot of \i\ (s = 1 
is fixed) and Q. m is presented in Fig. [6] With CMASS DR9 mea- 




0.26 0.28 Q 0.30 0.32 0.34 

m 

Figure 6. Confidence levels (lcr and 2cr) for joint fits to parameters fii (s = 
1) and il m . Solid lines show results for CMASS DR9 estimate of growth 
rate only, while dashed lines show results with previous estimates of growth 
rate added. In both cases the growth rate measurements are combined with 
CMB and SNela data. The dotted line shows the expected value in GR. 

parameter CMASS DR9 only CMASS DR9 + other f<rg 
y 0.75 ± 0.09 0.64 ± 0.05 

HX -0.24 ±0.11 -0.17 ±0.07 

Table 4. Constraints on parameters describing deviation from GR from 
CMB, SNela and growth rate data. 



surements we obtain /Ji = -0.24 ±0.11 when Q. m and H are 
marginalised over; after adding previous measurements of growth 
rate the constraint improves to = -0. 17±0.07. The GR predicted 
value of jii = 0.0 is again about 2<r from the best fit. 

A summary of constraints on MG parameters is presented in 
Table E] 

The constraining power of the growth rate data is strongly 
augmented by the fact that the initial amplitude of fluctuations at 
z ~ 1000 is very tightly constrained by CMB data. This CMB prior 
on the initial amplitude of the density fluctuations, which is degen- 
erate with the optical depth to the last scattering surface, allows 
us to convert our low redshift measurements of /(z)cr 8 (z) into pure 
measurements of growth rate f(z) that has a much stronger dis- 
criminative power. The detection of lensing effect on the power- 
spectrum of temperature and polarisation anisotropies from ongo- 
ing Planck satellite will enable a more robus t determination of the 
amplitude of primordial fluctuations (see e.g..lZaldarriaga & Sell aid 
Il998t IStompor & Efstafhioulll999l; lAmblard. Vale & Whitell2004 
which in turn will further boost the constraining power of RSD 
measurements on MG parameters. 

One of the key assumptions in the Reid et all d2012h analysis is 
that physics at recombination is very close to the "standard" model 
prediction, so that the CMB prior on the shape and amplitude of 
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Figure 7. Confidence levels (lcr and 2<x) for joint fits to parameters /i s and 
s. Grey areas show parts of the parameter space for which the effective 
Newton's constant at last scattering surface is modified by greater than 10 
per cent. 

the primordial fluctuations can be used. This assumption breaks 
down for the /j s parameterisation, for low values of .? that result in 
the rescaling of Newton's constant at earlier times. Modification 
of gravity at the last scattering surface will change the likelihood 
of observed CMB anisotropies and will not be consistent with the 
assumptions under which the measurements of iReid et al.l d2012l) 
were derived. 

Fig. UJ shows our constraints in /is space. Current CMB data 
provides a 10 per cent constraint on Newton's constant at high red- 
shift, which is ex pected to be red uced to a 1.5 per cent constraint 
with Planck data dGalli et al]|2009t) . The region of parameter space 
excluded by this limit is shown by the grey shading in Fig. UJ For 
the remaining area, we ignore the effect of the < 10 per cent change 
in Newton's constant on the primordia l power spectrum . 

Another important assumption in lReid et al. I d2012h is that the 
growth is scale independent. In many MG theories such as /(/?), the 
growth o f density fluctuation s is in fact strongly scale dependent 
(see e.g., Clifton etalj|2012h . Interpreting the growth rate data in 
the light of these models is not trivial since the measure ments were 
made assuming a scale independent growth (even though lReid et al.l 
120121 accounts for the scale dependence of power-spectrum due 
to nonlinear evolution). Testing these models with the growth rate 
data would require refitting the model to the measured anisotropic 
clustering while properly accounting for all the scale dependencies, 
which is outside of the scope of this paper. Our results could be con- 
sidered as constraints at a n effective scale, corresponding to those 
fitted bv lReid et al] fe012h . 

5.2 Constraining deviations from A 

For DE we consider parameters wo and w a , as described in Sec- 
tion |4.2| We will define the following models: 



WMAP + D v /r s 


+ anisotropic . 


\NMAP+D v /r s 




dlllbULlUpiC 
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Figure 10. Constraints from Rei d et all | |2012|) on the anisotropic galaxy 
clustering parameters fcrg and Fap at z = 0.57 are shown in red. The blue 
contour projects the combined constraints from WMAP and the spherically- 
averaged BAO position from BOSS, Dy/r, at z = 0.57; dashed green lines 
indicate the mean value of w in the region allowed by WMAP+Dy/r,; 
the black contour combines the anistropic measurement with the CMB and 
spherically-averaged BOSS BAO position, resulting in a factor of ~ 4 im- 
provement in the constraint on w. 



ACDM, where w = -1.0, w a = 0.0 and £\ = 0.0; 
wCDM, where wo is free but vv a = 0.0 and fi k = 0.0; 
vy w a CDM, in which both w and vv a are free and Q k = 0.0; 
OACDM, as ACDM but with free spatial curvature; 
OwCDM, as wCDM but with free spatial curvature. 

To demonstrate the constraining power and degeneracy directions 
resulting from our anisotropic measurements we first combine 
them with a WMAP7 prior on the shape of the primordial power- 
spectrum and amplitude only. We do this by using WMAP7 data 
to fit for a correlated joint likelihood of Q. m h 2 and the amplitude of 
matter fluctuations cr 8 at the CMB last scattering surface, with other 
parameters marginalized over. As we keep the amplitude of fluctu- 
ations fixed at the last scattering surface, the WMAP7 prior results 
in different values of <t s (z = 0) for different sets of parameters £l m 
and wo. Fig. [8]shows resulting constraints for ACDM and wCDM 
models. For the OACDM model we find n m = 0.308 ± 0.022 and 
100f2 k = 5.9 ± 4.8. In wCDM we obtain w = -0.91 ± 0.12. Fig.[8] 
clearly demonstrates how strongly complementary the RSD mea- 
surement of growth rate are to the geometric measurements. With 
our geometric and growth measurements alone, supplemented by 
WMAP7 shape and amplitude prior, we are able to constrain £l m 
with a 7 per cent precision and wo with a 14 per cent precision. 

Next, we combine BOSS measurements with the full WMAP7 
likelihood, so that the angular diameter distance to the last scatter- 
ing surface is also used. To strengthen these co nstraints we will als o 
combi ne the CM ASS DR 9 data with SNela dConlev et al]|201 ill , 
and H iRiess et alJ d201lh data sets. Table[5]presents the mean val- 
ues and la errors of cosmological parameters with different model 
assumptions and combination of these data sets. Fig.|9]presents the 
results for OACDM and wCDM models when CMASS DR9 re- 
sults are combined with full WMAP7 likelihood as well as either 
SNela or H data. 

When we assume a cosmological constant and fix the spatial 
curvature to zero (ACDM), we obtain a 5 per cent constraint on 
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Figure 8. Confidence levels (lcr and lcr) resulting from CMASS DR9 measurements and the WMAP7 shape and primordial amplitude prior. The solid lines 
show joint constraints from growth rate fcr% and geometrical measurements D\ and F, while the dotted lines show constraints from geometrical measurements 
only. The dashed line on left panel denotes the spatially-flat Universe, while the same line on the right panel denotes a cosmological constant (i.e. both lines 
show the locus of ACDM models). In both panels, w a = 0. 




0.24 0.26 0.28 0.30 0.32 0.34 0.24 0.26 0.28 0.30 0.32 0.34 0.36 

n n 

m m 



Figure 9. Confidence levels (lcr and 2<x) resulting from CMASS DR9 measurements combined with full WMAP7 likelihood, SNLS3 data of SNela and the 
estimate of Hubble constant. The solid line on the left panel denotes the spatially-flat Universe, while the same line on the right panel denotes a cosmological 
constant (i.e. both lines show the locus of ACDM models). In both panels, w a = 0. 
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Cosmological model 


Data set 




lOOOOi 


Ho 




w a 


ACDM 


CMB + CMASS + SNela 


0.285 ± 0.014 





68.9 ± 1.1 


.1 





ACDM 


CMB + CMASS 


0.291 ± 0.014 





68.5 ± 1.2 


.1 





ACDM 


CMB + CMASS + H 


0.281 ± 0.013 





69.5 ± 1.1 


.1 





OACDM 


CMB + CMASS + SNela 


0.281 ±0.014 


-9.2 ±5.0 


67.7 ± 1.3 


.1 





OACDM 


CMB + CMASS 


0.288 ± 0.017 


-8.5 ± -5.4 


67.4 ± 1.3 


.1 





OACDM 


CMB + CMASS + H 


0.277 ± 0.014 


-6.0 ±4.9 


68.8 ± 1.3 


.1 





wCDM 


CMB + CMASS + SNela 


0.292 ± 0.015 





68.0 ± 1.4 


-0.94 ± 0.05 





wCDM 


CMB + CMASS 


0.313 ±0.017 





65.9 ± 1.5 


-0.87 ± 0.05 





wCDM 


CMB + CMASS + H 


0.291 ± 0.015 





68.2 ± 1.4 


-0.93 ± 0.05 





OwCDM 


CMB + CMASS + SNela 


0.285 ± 0.017 


-8.2 ±5.5 


67.4 ± 1.5 


-0.98 ± 0.05 





OwCDM 


CMB + CMASS 


0.307 ± 0.022 


-3.9 ±6.8 


65.9 ± 1.6 


-0.90 ± 0.07 





OwCDM 


CMB + CMASS + H 


0.285 ± 0.018 


-3.7 ±5.7 


68.2 ± 1.5 


-0.95 ± 0.07 





Vt'oVl-aCDM 


CMB + CMASS + SNela 


0.280 ± 0.018 





68.8 ± 1.6 


-1.13 ± 0.12 


0.65 ± 0.36 


woWjiCDM 


CMB + CMASS 


0.313 ±0.037 





66.2 ± 2.8 


-0.86 ± 0.34 


-0.14 ± 1.04 


M>0W a CDM 


CMB + CMASS + H 


0.261 ± 0.037 





71.2 ±2.3 


-1.29 ±0.19 


1.02 ±0.48 



Table 5. Estimates of basic cosmological parameters for different models and combinations of data sets. 



Q m and a 2 per cent constraint on Hubble constant from the joint 
fit to CMB data, CMASS geometric measurements and SNela data. 
When the values of w a and Q.^ are fixed to zero, wo can be measured 
with a precision of 5 per cent. 



parameter 



most likely value lcr confidence level 



a 


0.16 


< 0.48 


Ml VM pl H 


1.0 


> 0.9 


<hlf 


0.11 


< 1.17 



These results ar e consistent with the results reported in 
lAnderson et al. who use only the position of the BAO peak 

and lSanchez et al] "i 20121) who use the full shape of the monopole. 
Our measurements are derived from the same sample but use 
additional information from the shape and the amplitude of the 
quadrupole. This allows us to get significantly stronger constraints 
with the CMASS data only (see Fig. [8}. When CMASS data is 
combined with other data sets this improvement is somewhat mit- 
igated but still present. When combined with CMB data and for 
the ACDM model our measurements improve constraints on Q. m 
by about 18 per cent relative to the BAO only results reported in 
lAnderson et al.l ( l2012l) . The central values are consistent within lcr. 
The biggest improvement is on the DE parameter w (— 0.87±— 0.05 
compared with -0.87 ± 0.24 in wCDM model), where we find a 
factor of four improvement over BAO only results even after com- 
bining the BAO measurements with WMAP7. Fig. 1 101 illustrates 
why including the information from anisotropic BOSS clustering 
improves the constraint on w so much, even though Fig. [4] shows 
that the marginalized one-dimensional constraints on /<x g and the 
Alcock-Paczynski parameter F are individually not very constrain- 
ing on w. The anisotropic BOSS constraint (red) is nearly per- 
fectly perpendicular to the degeneracy direction opened by w in the 
f<Ts(z = 0.57) - Fap(z = 0.57) plane for the data combination of 
WMAP7 with the BOSS measurement of D v (z = 0.57)/r s . When 
combined with CMB data our measurements of growth and geom- 
etry prefer wo > - 1 at about 2cr confidence level. When SNela or 
H data is added the agreement with ACDM becomes better. 



Table 6. Most likely values and lcr confidence level constraints on param- 
eters describing scalar field dark energy from our geometric and growth 
measurements and CMB shape and distance prior. 



5.3 Constraining the effective potential of scalar-field dark 
energy 

We also constrain parameters of scalar field dark energy models, as 
described in Section l4~3l for spatially-flat cosmological models. 

Fig. QT| shows constraints on the parameter a describing the 
potential of the slowly rolling scalar field and parameters M and 
/ describing the potential of PNGB field from our geometric and 
growth measurements combined with WMAP 7 shape and ampli- 
tude prior. For both models we assume a spatially-flat Universe. 
We also assume a uniform theoretical prior < // ^jM v \ < 1 and 
Ml -JM^Ho, where M v \ is the Planck mass. 

The resulting constraints on basic scalar field parameters are 
presented in Table [6] In both cases the ACDM limit of the scalar 
field theories provides a good fit to the data. The a parameter is 
constrained to be less than 0.48 at lcr confidence level, which is 
slightly better co mpared to the previously reported best constraint 
(Samushia 2009). M is needs to be at least 88 per cent of the Planck 
mass. For comparison previous studies used SNela data to constrai n 
M to be more than third of the Plank mass dDutta & Sorboll2008h . 
The ACDM limits of PNGB model (#,// = 0.0) and power-law 
model (a = 0) provide a good fit to data. 
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Figure 11. Confidence levels (lcr and 2<x) for joint fits to parameters fi„, and a for the inverse-power law model (left panel) and parameters M and / for 
the PNGB model (right panel) adopting the CMASS DR9 measurements of D\ and F and the WMAP7 shape prior. On the right panel, the shaded area 
corresponds to the region allowed at 2<x. Without adding the growth rate information the PNGB constraints are weak. The dotted line shows lcr contour from 
geometric measurements only, almost all of the remaining phase space is allowed at 2tr. 



5.4 Constraining model-independent properties of DE 

To reconstruct the model-independent properties of the late-time 
Universe, as des cribed in Section 14.41 we will use the geometric 
Rcid ct al. I d2012h measurements in combination with the previous 
AP constraints presented in Table [3] These estimates were derived 
using different fitting methods, measurements and range of scales. 
We ignore these differences in our analysis and will treat the re- 
ported measured values and their errorbars as Bayesian likelihoods 
of Dm and F at the redshifts of interest. 

If we assume that the Universe is well described by a FRW 
metric, the Hubble rate and expansion rate and the relationship 
with distances can be written without a reference to a particular DE 
model; they can be expressed purely as functions of a(z) through 

D A (z) = ^ f f. - , (22) 
Jo a(zf) 

(23) 



H(z) = a(z)d + z). 



where the scale factor a is normalised such that a = 1 today and 
dot denotes a time derivative. 

Different methods for the model-independent reconstruc- 
tion of the expansion history and p roperties of DE have been 
propo se d in the past (see e.g., lAlam. Sahni & Starobinsk 



20071: iTurner & Hutereil 120071 



& Diorgovski 2007 



insky 
20071 



Shafieloo & Clarksonl l20ld ; Crittenden et al .1 I2012L and refer- 



ences therein). The most popular approach is to parametrize the 
dimensionless coordinate distance in the integral in Eq. (I22t by a 
polynomial and express the sca le factor through its first and second 
deriva tives. iBlake et al. Il l2011ch used the method of Shafieloo et al. 
J2006I) to reconstruct the comoving distance from SNela data, 
which, combined with their measurements of F, enabled an 
estimate of H(z). 

We will adopt an alternate approach that does not require tak- 



ing derivatives of the quantities reconstructed from data, which has 
been shown to introduce artificial oscill ations in reconstructed vari- 
ables dLazkoz. Salzano & Sendral2 012). We assume that a 3 H 2 (a) = 
cut 2 = y(a) is a smooth function of the scale factor, so that it can be 
approximated by a third order polynomial. We have checked that 
this approximation holds to a percent level for all conventional DE 
models for wide range of redshifts; for ACDM this expansion is 
exact up to the radiation dominated era. Eqs. ((23} can be rewritten 
in terms of y(a) as 



D A (a) 



= ca I 



a' 3/2 da' 



H(a) = 




(24) 
(25) 



We map y(a) = ^ y\d into D\ and F using Eqs d24t - i25l 

;=o 

and compare the results to geometric measurements. We compute a 
likelihood of polynomial coefficients of the expansion of -£(a,) and 
map it onto the likelihood of deceleration parameter q{z) which is 
defined in Eq. dl9l > by 



£\q) = £(y[a,]). 



(26) 



The left panel of Fig.[T2]presents the results of the reconstruc- 
tion of q(z) using our geometric measurements from the CMASS 
DR9 data combined with previous similar geometric measurements 
from 6dFGRS, SDSS-II and WiggleZ surveys. The reconstruction 
uses only the geometric measurements, shape prior from WMAP7 
data and theoretical priors of spatially-flat FRW metric and smooth- 
ness of a 3 H 2 (a). The recovered deceleration parameter is less than 
zero at redshifts z < 0.3 and more than 3cr from q = 0.5 as pre- 
dicted by EdS model. La mpeitl et alj d201fj|) were able to derive a 
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Figure 12. Confidence levels (ltr and 2<x) for the deceleration parameter as a function of redshift and Om(z) reconstructed from the compilation of geometric 
measurements in tables [2] and [5] Ho is marginalised over with an HST prior. The dotted line in the left panel demarcates accelerating expansion (below the 
line) from decelerated expansion (above the line). The dashed line in both panels shows the expectation for an EdS model. 



much stronger constraint of q = -0.34 ± 0.18 using SNela data 
only, however they had to assume constant q while we allow q to 
have non-zero higher derivatives. 

The right panel of Fig. [represents a similar reconstruction of 
Om(z) parameter. The reconstructed Om(z) is consistent with be- 
ing constant and is more than 5<x from the EdS predicted value 
of Om(z) = 1. Our reconstr ucted q(z) and Om( z) are consistent 
with those reconstructed by iBlake et all d201ld) : our results are 
smoother as a function of redshift because of the different recon- 
struction method used. 

Our geometric measurements can be used to derive an estimate 
of Hubble expansion rate at z = 0.57, since 



5.5 Is there an evidence for deviations from flat LCDMGR? 



H(z) 



2U/3 



Dv 



(27) 



This measurement of H(z) combined with an independent measure- 
ment of H can be used to estimate the increase in the fractional 
energy density of DE with respect to nonrelativistic matter from 
z = 0.57 to present day. If we assume a flat FRW background and 
two dominant components fl m and Q DE then 



p DE (z) = H(z) 2 /H 2 -£l m g+z? 
p m (z) fi m (l+z) 3 



(28) 



Comb ining our estimate of H(z = 0.57) with the iRiess et all 
d201 lh measurement of H and a WMAP 7 prior of £2 m /i 2 = 
0.1334 ± 0.0056 we find p DE (z = 0)/p m (z = 0) = 3.09 ± 0.32 and 
p DE (z = 0.57)/p m (z = 0.57) = 0.67±0.13. This means that the frac- 
tional contribution of DE to the total has increased by a factor of 
~ 4.7 since z = 0.57, lending support to the notion that expansion 
dominated by DE is a relatively recent occurrence. 



As discussed in iReid et alj J2012h the observed galaxy clustering 
when combined with CMB data is fully consistent with LCDMGR 
expectations. Nevertheless, several of our fits to one-parameter ex- 
tensions of flat LCDMGR (y, fix, and w ) indicate a ~ 2cr prefer- 
ence for values away from the fiducial. The data, however, demand 
these extra parameters at a much lower significance than 2cr: the 
difference in the best fit^ 2 values between flat LCDMGR and these 
one-parameter model extensions is only ~2.5. A similar effect has 
been observed for the effective number of Neutrinos (N e g = 3 in the 
standard model of particle physics), where the Bayesian posterior 
likelihood resulting from the CMB data has been shown to pre- 
fer jVe ff > 3 at high confidence level. The lGonzalez-Morales et all 
J201 ]]), however, showed that in this case the priors imposed on 
cosmological parameters have a strong effect on the posterior like- 
lihood of A'ejf and after removing this prior-dependency the prefer- 
ence for the deviation from N s g = 3 is much lower. 



6 CONCLUSIONS 



We have used the IReid et al.1 J2012h measurements of angular 
distance, Hubble expansion rate and growth rate derived from 
the anisotropic clustering of BOSS CMASS DR9 galaxies to 
place constraints on deviations from the standard cosmological 
model that assumes a ACDM background with structure formation 
driven by GR. The geometric measurements of D v and F are 
compl ementary to similar m easurements from the BAO peak po- 
sition (Anderson et al.l 20121) and the full shape of the correlation 
function dSanchez et alj|2012r) and strengthen existing constraints 
on parameters describing the time-dependence of DE energy 
density. The RSD measurement of fa% was shown to provide an 
additional constraint on the parameters describing deviations from 
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GR and helped to significantly tighten DE constraints derived 
from geometric measurements. We now highlight our findings 
by using them to answer fundamental questions about our Universe. 

How much do RSD measurements enhance the geometric measure- 
ments ? 

When GR is assumed the RSD measurements of growth break 
parameter degeneracies present when using purely geometric 
measurements, and consequently significantly tighten constraints 
on basic cosmological parameters. The addition of growth rate 
information improves constraints on £l m by 18 per cent relative 
to the case where only BAO peak position data is used. The 
improvement is most dramatic for the w parameter constraints, 
which improve by a factor of four compared to BAO only results 
(see Table [5} . 

Does GR provide a good description of data? 
When ACDM is assumed our measurements of growth and 
geometry show a 2<x preference for models in which the growth 
of structure is weaker compared to GR. Adding previous RSD 
measurements at other redshifts brings the best-fit closer to GR 
value but still prefers somewhat weaker growth. This results from 
the fact that most RSD measurements of f(z)crg(z) with high 
signal-to-noise are below GR predictions (see Fig.[T}. 

Does the Cosmological Constant provide a good description of 
data? 

Assuming GR, our measurements of growth and geometry show 
a 2cr preference for w > -1. When combined with SNela data, 
the best-fit is closer to the cosmological constant model and the 
discrepancy is lowered to about lcr. It should be noted that, for DE 
as well as GR, the 2a preference means a preference in terms of 
relative Bayesian likelihood. 

How well can the DE scalar field potential be constrained? 
We demonstrated that our measurements of growth and geometry, 
when combined with CMB information, provide strong constraints 
on scalar field DE model parameters. The constraints obtained are 
better than previously reported from different combinations of data 
sets. The flat potential (Cosmological Constant) provides a good 
fit to data. 

Has the expansion of the Universe accelerated recently? 
We employed a model-independent approach that relies on very 
few general assumptions to reconstruct the "deceleration param- 
eter" at low redshifts. We showed that current AP measurements 
provide 2 to 3cr evidence for the accelerated Universe at low red- 
shifts. Much stronger measurements of q(z) exist in the literature, 
but they depend on additional assumptions about the nature of DE. 

Did DE emerge as a dominant component only very recently? 
Our estimate of H(z = 0.57) combined with an estimate of Ho 
suggest that energy density of DE relative to nonrelativistic matter 
was about 4.5 times lower at z = 0.57. 

Is the standard ACDMGR model still valid? 
Measurements of growth and geometry from CMASS DR9 sample 
allow for a very strong test of MG and DE. When combined with 
WMAP7 data they show a 2cr preference for either weaker grav- 
ity or w > -1. However, in terms of absolute ,^ 2 /dof, the simple 
ACDMGR model still provides a good fit to the data. 
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APPENDIX A: NONGAUSSIANITY IN POSTERIOR 
LIKELIHOOD 

Figure lATl shows the marginalized one-dimensional posterior like- 
lihood of parameters fcrg, D A and H used in our analysis. Our mul- 
tivariate Gaussian approximation to the likelihood works very well 
everywhere except for the tails of the likelihood where it deviates 
from the empirical likelihood by not accounting from for small pos- 
itive skewness in /erg and H. We have checked that the derived 
likelihood of modified gravity and dark energy parameters does not 
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Figure Al. Marginalized one-dimensional posterior likelihood of parameter /o"8, Da and H from Reid et al. 12012) used in our analysis (The likelihood is in 
arbitrary units). The bars show the real distribution, while the solid line represents the Gaussian approximation that we use. 



change at 1 and 2cr confidence levels when the Gaussian approxi- 
mation is used instead of the whole likelihood. 
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